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Fig. 1. Mapped outlines of slump and debns avalanche aprons. outline in, he bottom right is that of a possible 
lava flow from the collapsed margin dome on Sapas Mons. 


On Earth, landslides on volcanic edifices can be triggered by a 
number of different processes, including those occurring as a result 
of aseismical crustal deformation, such as oversteepen.ng of slopes 
due to deformation (possibly resulting from dyke emplacement o 
magma rise), overloading of the slope (by lavas), excess weight at 
the top of the slope (due to a large cone or a large area of summit 
lava), removal of support by explosions on the flanks, and caldera 
collapse. Failure occurring coseismically can result from structural 
alteration of the constituent parts of the slope leading to failure, 
dislodgement of otherwise stable slopes , and fault movement resum- 
ing in an increased slope angle [4], Seismic pumping may also be a 
major control on slope stability during an earthquake [5]. 

On Venus, similar processes may operate. The high ambient 
temperatures may result in development of a weak carapace, which 
in turn may allow relatively rapid dome growth to occur. If the 
effusion rates are high, as suggested by the size of the features, then 
oversteepening would be a likely consequence resulting ,n failure 
and collapse. Landslide scars may be modified by continued dome 
growth The existence of fractures around the base of some of the 
collapsed domes and of debris aprons cut by fractures suggests that 


rc has been seismic activity and surface deformation occurring 
ing the period of modification of the dome. 
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CED-VALENCE IRON MINERALS ON VENUS: F^-FV* 
[DES AND OXY-SILICATES FORMED BY SURFACE- 

VIOSPHEREINTERACTIONS. Roger G. Bums and D'Arcy 
Straub. Department of Earth. Atmospheric and Planetary Set- 
cs, Massachusetts Institute of Technology, Cambridge MA 


Background: The oxidation state and mineralogy of iron on the 
hot surface of Venus are poorly understood [1-3], despite qualita- 
tive in situ measurements of oxygen fugacity during the Venera 13/ 
14 missions [4), some reflectance spectral data derived from the 
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rable temperatures, the crystallization of babingtonite requires 
more hydrous conditions , lower C0 2 , and slightly higher 0 2 fugact- 
ties in the fluid phase than ilvaite. Since similar temperatures, C0 2 
pressures , and oxygen fugacities induced within skam deposits exist 
on Venus, ilvaite and perhaps babingtonite could have also formed 
on the surface of this planet by the interaction of the venusian 
atmosphere with extruded basaltic rocks. One factor that might 
mitigate against the formation of these calcic Fe*--Fe* silicates on 
Venus however, are the high abundances of Mg and A1 measured 
during the Venera 13/14 [32] and Vega 2 [33] missions. The Mg 
and Al 1 * cations are not accepted into the crystal structures of ilvaite 
and babingtonite. 

Discussion: Although magnetite is generally regarded to be 

the predominant ferric-bearing mineral on Venus, other muted- 
valence Fe 2 *-Fe 5+ minerals known to exist on the surface of fearth 
could be stable in the venusian atmosphere. Thus, in addition to 
laihunite(whichis probably metastable) and ilv aite and babingtonite 
(both of which may be found in rocks depleted of Mg and Al), oxy- 
amphiboles and oxy-micas may also be major constituents of the 
venusian surface. The opacities and high electrical conductivities of 
such mixed-valence Fe^-Fe 3 * silicate minerals, the properties of 
which resemble magnetite [34], may also contribute to high radar- 
reflectivity regions in the highlands of Venus [35]. 
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THERMAL BUOYANCY ON VENUS: PRELIMINARY RE- 
SULTS OF FINITE ELEMENT MODELING. J. D. Burt and 
J. W. Head, Department of Geological Sciences, Brown University, 
Providence RI 02912, USA. 


Introduction: Enhanced surface temperatures and a thinner 
lithosphere on Venus relative to Earth have been cited as leading to 

increased lithospheric buoyancy. This would limit[l] or prevent [2] 

subduction on Venus and favor the construction of thickened crust 
through underthmsting. Underthrusting may contribute to the for- 
mation of a number of features on Venus. For example, Freyja 
Montes, a linear mountain belt in the northern hemisphere of Venus , 
has been interpreted to be an orogenic belt [3] and a zone of 
convergence and underthmsting of the north polar plains beneath 
Ishtar Terra, with consequent crustal thickening (4] . Such mountain 
belts lie adjacent to regions of tessera and contam evidence of 
volcanic activity. Tessera is also considered to consist of thickened 
crust [5] and crustal underthmsting is one possible mode for its 
formation [4]. Models for the formation of the mountain belts and 
associated features must then explain compressional deformation 
and crustal thickening, as well as melt production. 

In order to evaluate the conditions distinguishing between 
underthmsting and subduction, we have modeled the thermal and 
buoyancy consequences of the subduction end member. This study 
considers the fate of a slab from the time it starts to subduct, but 
bypasses the question of subduction initiation. Thermal changes in 
slabs subducting into a mantle having a range of initial geotherms 
are used to predict density changes and thus their overall buoyancy. 
Finite element modeling is then applied in a first approximation of 
the assessment of the relative rates of subduction as compared to the 
buoyant rise of the slab through a viscous mantle. 

Subduction Model: In the model, slabs, having a thickness set 

by 90% of the basalt solidus, subduct at a 45° angle into the mantle. 
The initial geothemns match surface thermal gradients of l0°C/km, 
l5°C/km, and 25°C/km [6]. Slabs heat via conduction, crustal 
radioactivity, phase changes, and adiabatic compression. Phase 
changes involving the conversion of basalt to eclogite al depths of 
60 to 1 60 km and then enstatite to forsierite plus stishovite between 
260 and 360 km generate 0. 1 3 x 1 (h 5 crgs/cm 3 s and 0 .3 6 x 1 (h 5 ergs/ 
cm 3 s respectively [7]. The slab radiogenic heat production is 2 .63 x 
10-7 ergs/g s [81. Adiabatic compression adds 0.5°C per kilometer 
of depth. Convergence rates ranged from 5 mm/yr to 100 mm/yr. 

The thermal evolution of the slab is followed using a finite 
difference technique [7,9]. The model region measures 800 km 
horizontally by 400 km deep. Processing ends when the slab tips 
reach a 300-km depth, imply ing time intervals of 1 0 m.y . to 1 00 m.y . 
Slab density changes derive from the thermal results through 
calculation of the thermal expansion = 3 x 10~ 5 /°K W) and the 
effects of pressure (b = 1 x 10- 3 /kbar [81) on initial densities set for 
zero pressure and temperature. The assumed initial density structure 
includes a 10-km or 25-km basaltic crust (density = 3.0 g/cm 3 ), a 
corresponding 25-km or 65-km-thick depleted mantle zone (den- 
sity _ 3.295 g/cm 3 ), and an underlying undepleted mantle (den- 
sity = 3.36 g/cm 3 ). Density changes due to the phase transitions are 
also included. Results take the form of density distributions within 
the model region. 

Finite element modeling is then employed to gauge the rale at 
which a slab having the density structure derived above will move 
through a viscous mantle. Buoyancy body forces are applied to a 
slab having a viscosity of 10 21 Pa s surrounded by a mantle with a 



